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ELASTIC DISPLACEMENTS FOR VARIOUS EDGE-CRACKED PLATE SPECIMENS 
by Bernard Gross, Ernest Roberts, Jr., and John E. Srawley 
Lewis Research Center 

SUMMARY 

The relative displacement per unit load of two conjugate points is used as a quantita- 
tive indicator of crack extension in plane strain fracture toughness Kj c tests. The nec- 
essary displacement data are presented here in dimensionless form for five types of 
single-edge-crack specimens: three-point bending, pure bending, remote axial tension, 
and eccentric tension of compact rectangular and tapered varieties. The results were 
obtained by a boundary collocation method of elastic analysis and are highly precise in 
themselves. The accuracy with which they apply to actual specimens depends on the ex- 
tent to which the assumed boundary conditions are equivalent to the actual load distribu- 
tions. Reasonably good agreement was obtained with three sets of existing experimental 
data, and this agreement is sufficient for Kj c test purposes. For gage points located on 
the specimen edge, the results show low sensitivity to small variations in gage location. 


INTRODUCTION 

In plane strain fracture toughness Kj c testing, the relative displacement per unit 
load Vy/P of two conjugate points is used as a quantitative indicator of crack extension 
during the test (refs . 1 to 3) The relation between Vy/P and relative crack length a/W 
depends on the type of specimen employed. The relations for various specimens are de- 
termined by linear elastic strain analysis, or by direct experimental measurements in 
which slots are used to simulate cracks. This report presents Vy/P data obtained by a 
boundary collocation method of analysis for several types of single -edge -crack plate 
specimens. The use of these specimens in Kj c testing has been explained by Brown and 
Srawley (refs . 1 and 2) . 

The following types of loading were investigated (figs. 1 and 2): three-point bending 
having a span to width ratio of 4, pure bending, remote axial tension, and eccentric ten- 




Three-point bend Pure bend specimen 

specimen (S/W = 4) 


(a) Three types of loading. 



(b) Common geometric configuration. 

Figure 1. - Bending- and axial-tension plate specimens. 
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sion (compact rectangular and tapered varieties of specimens) . The results are given in 
the form of tables of a dimensionless displacement coefficient as a function of the major 
variable, relative crack length a/W, and of the secondary variables x/W and y/W 
which represent the gage point locations in relation to the crack tip. 

Displacement measurement for detection of crack extension in a Kj c test should not 
be confused with measurements of specimen compliance (reciprocal stiffness) for deter- 
mination of energy release rates (ref. 1). The intent of a compliance experiment is to 
determine the work done by the loading forces, and the displacement measured must be 
chosen appropriately. For evaluation of crack extension, any convenient gage points can 
be used for the displacement measurement. The most sensitive positions are those close 
to the crack, which are precisely the least suitable for compliance when the specimen is 
remotely loaded. However, for the eccentric tension specimens discussed in this report, 
the gage point positions that are convenient for crack extension detection are also suitable 
for compliance determinations. 
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Figure 2. - Eccentric tension plate specimens. 


SYMBOLS 


a 

B 

d 2n ,d 2n-l 

E 

e 

H 

h p 

K 

K lc 

M 

P 


crack length 
specimen thickness 

coefficients of Williams stress function 
Young’s modulus 

distance from wedge tip to line of load application 

uniform depth of nontapered split arm 

depth of tapered split arm at load line 

stress intensity factor of crack tip elastic stress field 

plane strain fracture toughness, measured in terms of the opening mode 

1 /2 

stress intensity factor Kj in units of (Stress) x (Length) ' (see refs. 

1 and 2) 

bending moment 
load 

polar coordinates referred to crack tip 
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S span 

U displacement in radial direction 

u x displacement in x-direction 

V displacement in tangential direction 

Vy displacement in y-direction between ±y locations 

W specimen width 

x, y Cartesian coordinate system referred to crack tip 
9 polar coordinate referred to crack tip 

v Poisson's ratio 

(p harmonic function in displacement equation 

X Airy stress function 


ANALYSIS 

The method of analysis is described in some detail by Gross and coauthors (refs. 

4 to 8). Its earliest use in solid mechanics is by Barta (ref. 9). Considerable detail is 
given by Green (ref. 10) and Howland and Knight (ref. 11). The method is called various- 
ly "boundary collocation" and "point matching". Briefly, it consists of truncating a 
series solution to the appropriate partial differential equation, and making use of the 
boundary values at a finite number of points to evaluate its coefficients. It can be shown 
(ref. 12) that the biharmonic equation in terms of an Airy stress function properly de- 
scribes the plane elastic problem. It is convenient to use the Williams stress function 
for our analysis (refs. 13 and 14). It is an Airy stress function, it identically satisfies 
the biharmonic equation, and it identically satisfies the boundary conditions along the 
crack surface. 

The equation to be solved is 


v 4 x = o 

where the origin of the polar coordinate system is the crack tip (fig. 3). The solution 
given by Williams is 
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x(r, 0) = 


OO 



n=l, 2, 3. 



l d n+1/2 
°2n-l r 


-cos 



2n - 3 
2n + 1 


cos 



+ ( -1 ) nd 2n rn+1 [" cos ( n 


1)0 + cos(n + 




Figure 3. - Geometry defining displacements in x- and y-directions. 


The displacements in terms of the stress function and a harmonic function cp(r, 9) are 
given by Williams (refs. 13 and 14) and Coker and Filon (ref. 15). In infinite series form 
they become 


OO 



2juV(r, 9) = - A + (1 - a)r 2 ®S? 

r 9r d9 


2juU(r, 9) = - & + (1 - a)r iS? 

ar 30 


5 



where for plane stress a = v/(l + v), for plain strain a = v, and for either plane stress 
or plane strain ji = E/[2(l + v)]. The resulting displacement equations are 


V(r, 0) = — ) 

2 n L^l 


'i-ll n d r n-l/2 
l) a 2n-l r 


n=l, 2, 3. . . 


- + n - 4a) sin (n - - ) 0 - (— - sin (n + - |0| 


(-l) n d 2 n r n j(n - a)sin(n - 1)0 - (n + l)sin(n + 1)0] 


U(r, 0) = — 

2p 


n=l, 2, 3. . . 


i I (■ l>nrn ’ 1/2dan - 1 16 - ■ 4 j cos f ■ + f - i) cos ( n + 1 ) 


+ ( _ l) n+ ^2 n [(3 - n - 4a)cos(n - 1)0 + (n + l)cos(n + 1) 0 ]^ 


For the special case of 0 = -n, it can be deduced from these equations that the dis- 
placements for plane stress are independent of Poisson’s ratio v, and that the displace- 
ments for plane strain are equal to (1 - v ) times those for plane stress. Since this 
simplification applies to nearly all the results obtained, it was convenient to compute the 
results for plane stress. 

The displacements in the x- and y-directions (fig. 3) are 

u = U cos 0 - V sin 0 

A 

Vy = 2(U sin 0 + V cos 0) 

It is to be noted that v y as defined in the preceding expression is the displacement in- 
dicated by a gage mounted across the crack, that is, the relative displacement of a pair 
of conjugate points. 

The number of terms in the series, and hence the number of boundary points satis- 
fied, was progressively increased until negligible changes occurred in the values of the 
displacement. The relation between calculated displacement and the number of terms in 
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Number of matrix coefficients used in computing v y BEW/6M 


Figure 4. - Typical variation of dimensionless displacement coefficient v y BEW/6M 
as function of number of coefficients used in computation. Displacements are 
for three-point bend specimens computed at specimen crack edge. 


the series is shown graphically in figure 4 for one specimen configuration. 

All quantities used in the analysis were made dimensionless. Hence, all crack 

lengths are relative to specimen width a/W, all coordinates are relative to specimen 

width x/W and y/W, and all displacements are relative to elastic modulus, specimen 

thickness, and either load or moment per unit depth W, v BE/P and v EBW/6M. 

y y 

RESULTS AND DISCUSSION 

Dimensionless displacement coefficients for plane stress are tabulated for values of 
the relative crack length a/W up to 0.7 in tables I to VI. The results for three-point 
bend specimens S/W = 4 (tables I and n) are more detailed than the others because these 
specimens are of particular current interest for standardized Kj c measurement (ref . 2). 
Tables I and II show the displacements for different gage point locations, respectively; 
along the crack, and along the specimen edge normal to the crack. The displacement is 
considerably more sensitive to position along the crack than along the specimen edge 
normal to the crack. For this reason, it is desirable in practice to locate the gage on 
the edge of the specimen. However, in Kj c measurement, the factor needed is the de- 
rivative of the logarithm of the displacement coefficient with respect to the logarithm of 
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the relative crack length (ref. I) 1 , and this factor is less sensitive to gage location than 
is the displacement coefficient itself. Therefore, the results obtained in K Ic tests will 
show low sensitivity to small variations in gage location. 

As a matter of interest, the displacements parallel to the crack u x are given in 
table I as well as the displacements normal to the crack. These lengthwise displacements 
are comparatively small and are not of any particular interest. 


TABLE I. - DIMENSIONLESS DISPLACEMENTS ALONG CRACK FOR 
THREE-POINT BEND SPECIMENS (S/W = 4) (PLANE STRESS) 


Gage location 



a/W 



x/W 

y/w 

0.3 

0 . 

5 

0.7 




EBv y W/6M 

EBu x W/6M 

EBv y W/6M 

EBu x W/6M 

EBv y W/6M 

EBu x W/6M 

£ 

CM 

00 

N 

d 

±0 

0.316 

0.0012 

0.722 

-0.0024 

1.948 

-0.0230 

-a/16W 



.450 

.0025 

1.041 

-.0045 

2.876 

-.0436 

-a/8W 



.644 

.0052 

1.525 

-.0080 

4.391 

-.0789 

-a/4W 



.933 

.0107 

2.298 

-.0120 

7.050 

-.1303 

-a/2W 



1.377 

.0218 

3.611 


12.019 

-.1802 

-a/W 



2.099 


5.972 


21.597 



The last line in table n lists experimental results by Fisher (ref. 3) for comparison. The 
experimental results are somewhat lower than the computed results for plane stress, but 
not as much as by the factor (1 - v ) for v = 0. 3. Thus, the experimental results are 
bracketed by the computed results for plane stress and those for plane strain (to a close 
approximation). This would be expected since the region near the crack tip in the actual 
specimen approaches a state of plain strain, whereas regions remote from the crack tip 
are in a state of plane stress. The two-dimensional analysis is not capable of producing 
closer agreement with the experimental results. Experimental results are also com- 
pared with computed results in table HI for pure bending, and the differences are similar 
to those for three -point bending. The displacements are greater for pure bending than 
for three -point bending because of the different bending moment distributions. 

Table IV gives results for single-edge-cracked specimens under remote axial ten- 
sion on the assumption of uniform stress distribution at a distance not less than 0. 8 W 
from the crack (ref. 4). These results are intended to apply- to pin-loaded specimens 
with pin centers not less than 3 W apart. 


*For example, for remote axial tension, the factor is [d log(VyEB/P)/d log(a/W)| 
which is equal to[(a/W)/(v EB/P)J[d(v EB/P)/d(a/W)j as in fig. 40 of ref. 1. 
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TABLE II. - DIMENSIONLESS DISPLACEMENTS AT EDGE FOR THREE- 


POINT BEND SPECIMENS (PLANE STRESS) WITH COMPARATIVE 


EXPERIMENTAL RESULTS (S/W = 4) 


Gage location 



a/W 



Source 

x/W 

y/w 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 






EBv y W/6M 




-a/W 

±0 

1.184 

2.091 

3.520 

5.973 

10.733 

21.597 

Collocation 


±.10 

1.217 

2.109 

3.529 

5.987 

10.751 

21.625 

Collocation 


±.20 

1.453 

2.293 

3.696 

6.154 

10.940 

21.865 

Collocation 


±.50 

2.385 

3.082 

4.388 

8.817 

11.632 

22.656 

Collocation 


±.10 


1.94 

3.46 

5.60 

9.74 

19.9 

Fisher, et al. 









(ref. 3) 


TABLE HI. - DIMENSIONLESS DISPLACEMENTS AT EDGE FOR PURE BEND 
SPECIMENS (PLANE STRESS) WITH COMPARATIVE 
EXPERIMENTAL RESULTS 


Gage location 



a/W 



Source 

x/W 

y/w 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 






EBVyW/6M 




-a/W 

O 

O -H 
-H 

1.260 

1.20 

2.241 

2.17 

3.763 

3.69 

8.368 

6.00 

11.354 
10. 80 

22. 575 
21.40 

Collocation 
Fisher, et al. 
(ref. 3) 


TABLE IV. - DIMENSIONLESS DISPLACEMENTS AT EDGE 
FOR REMOTE AXIAL TENSION SPECIMENS 
(PLANE STRESS) 


Gage location 

a/W 

x/W 

y/w 

0.2 0.3 0.4 0.5 

0.6 

0.7 



EBVy/P 

-a/W 

±0 

1.440 2.806 5.217 9.881 

19.900 

44.100 


I I 


III Hill 



Tables V and VI give results for various compact tension specimens which are ec- 
centrically loaded at positions close to the crackline and to the specimen edge (ref. 8). 
Table V includes a comparison with experimental results from Bush and Wilson (ref. 16), 
which are in satisfactory agreement with the computed results in spite of the consider- 
able difference of gage point location. It would be expected, however, that this type of 
specimen would be particularly insensitive to gage location along the specimen edge. 


TABLE V. - DIMENSIONLESS DISPLACEMENTS AT LOAD LINE 
FOR ECCENTRICALLY LOADED COMPACT RECTANGULAR 
TENSION SPECIMENS (PLANE STRESS) 


Gage location 

W/H 

a/W 

Source 

x/W 

y/w 

0.3 

0.4 

0.5 

0.6 

0.7 

EBVy/P 

-a/W 

±0 

5/4 

11.813 

19. 277 

32. 140 

57. 335 

116. 031 




5/3 

14. 266 

22. 968 

37.015 

63. 320 

122. 050 




20/11 

15. 125 

25.058 

39. 877 

66. 926 

126. 699 




25/12 

18.075 

29.473 

46. 229 

75. 234 

136. 966 


x/W 

y/w 

W/H 



a/W 



Source 




0.333 

0. 389 

0. 444 

0. 500 

0. 556 


-a/W 

;±o 

9/4 

23. 392 

30. 480 

40.057 

51. 145 

65. 964 

Collocation 


±. 875 

9/4 

24.96 

30.96 

38. 52 

49. 20 

63.36 

Bush, Wilson 









(ref. 16) 


TABLE VI. - DIMENSIONLESS DISPLACEMENTS AT LOAD 
LINE FOR ECCENTRICALLY LOADED TAPERED 
TENSION SPECIMENS (PLANE STRESS) 


Gage location 

H p/e 

W/e 

a/W 

x/W 

y/w 

0.2 

0.3 

0.4 

0. 5 - 

0.6 

0. 7 

EBVy/P 

-a/W 

±0 

0.2 

3.2 

162.8 

312. 7 

474. 2 

641.9 

810.0 

985.4 



.3 

4.2 

97. 3 

170.8 

242.4 

313. 6 

387. 5 

481.0 



.4 

5.0 

67.0 

107.8 

148.3 

187. 3 

234.9 

306.3 
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CONCLUSIONS 


Boundary collocation is a satisfactory procedure for computing elastic displacements 
(per unit load) for single-edge-crack specimens. The results are highly precise and in- 
expensive of computer time, so that variations of specimen shape and gage location can 
be explored at little cost. To do the same experimentally would be quite expensive. The 
accuracy with which the results apply to actual specimens depends, of course, on the ex- 
tent to which the assumed boundary conditions are equivalent to the actual load distribu- 
tions. 

Where the present results could be compared with existing experimental results 
shown in the tables, the agreement is reasonably good and is adequate for plane strain 
fracture toughness Kj c test purposes. The results confirm the expected low sensitivity 
to small variations in gage point location when these points are on the specimen edge. 

The sensitivity to change of position of gage point location along the crack is greater. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 13, 1967, 

731-21-03-01-22. 


REFERENCES 

1. Brown, William F. , Jr. ; and Srawley, John E. : Plane Strain Crack Toughness 

Testing of High Strength Metallic Materials. Spec. Tech. Pub. No. 410, ASTM, 

1966. 

2. Srawley, John E. ; Jones, Melvin H. ; and Brown, William F. , Jr. : Determination of 

Plane Strain Fracture Toughness. Mat. Res. & Standards, vol. 7, no. 6, June 

1967, pp. 262 - 266. 

3. Fisher, Douglas M. ; Bubsey, Raymond T. ; and Srawley, John E. : Design and Use 

of Displacement Gage for Crack- Extension Measurements. NASA TN D-3724, 1966. 

4. Gross, Bernard; Srawley, John E. ; and Brown, William F. , Jr. : Stress-Intensity 

Factors for a Single -Edge-Notch Tension Specimen by Boundary Collocation of a 
Stress Function. NASA TN D-2395, 1964. 

5. Gross, Bernard; and Srawley, John E. : Stress -Intensity Factors for Single-Edge- 

Notch Specimens in Bending or Combined Bending and Tension by Boundary Colloca- 
tion of a Stress Function. NASA TN D-2603, 1965. 


11 



6. Gross, Bernard; and Srawley, John E. : Stress-Intensity Factors for Three-Point 

Bend Specimens by Boundary Collocation. NASA TN D-3092, 1965. 

7. Gross, Bernard; and Srawley, John E. : Stress -Intensity Factors by Boundary Collo- 

cation for Single-Edge-Notch Specimens Subject to Splitting Forces. NASA TN 
D-3295, 1966. 

8. Srawley, John E. ; and Gross, Bernard: Stress Intensity Factors for Crackline- 

Loaded Edge-Crack Specimens. Mat. Res. & Standards, vol. 7, no. 4, Apr. 1967, 

pp. 155-162. See NASA TN D- 3820, 1967. 

•• •• •« 

9. Barta, J. ; Uber die Naherungsweise Losung Einiger Zweidimensionaler Elasti- 

zitatsaufgaben. Z. Angew. Math. Mech. , vol. 17, no. 3, June 1937, pp. 184-185. 

10. Green, A. E. : General Bi -Harmonic Analysis for a Plate Containing Circular Holes. 

Proc. Roy. Soc. London, ser. A, vol. 176, no. 964, Aug. 28, 1940, pp. 121-139. 

11. Howland, R. C. J. ; and Knight, R. C. : Stress Functions for a Plate Containing 

Groups of Circular Holes. Phil. Trans. Roy. Soc. London, ser. A, vol. 238, 

1940, pp. 357-392. 

12. Sokolnikoff, Ivan S. : Mathematical Theory of Elasticity. Second ed. , McGraw-Hill 

Book Co. , Inc. , 1956. 

13. Williams, M. L. : On the Stress Distribution at the Base of a Stationary Crack. J. 

Appl. Mech., vol. 24, no. 1, Mar. 1957, pp. 109-114. 

14. Williams, M. L. : The Bending Stress Distribution at the Base of a Stationary Crack. 

J. Appl. Mech. , vol. 28, no. 1, Mar. 1961, pp. 78-82. 

15. Coker, E. G. ; and Filon, L. N. G. : A Treatise on Photo -elasticity. Second ed. , 

Cambridge University Press, 1957. 

16. Bush, A. J. ; and Wilson, W. K. : Determination of Energy Release Rate for Biaxial 

Brittle Fracture Specimen. Rep. No. WERL-8844-2, Westinghouse Electric Corp. , 
Aug. 1964. (Available from DDC as AD-611659. ) 


12 


NASA- Langley, 1967 32 E-4016 



"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof ” 

—National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and technical information considered 
important, complete, and a lasting contribution to existing knowledge. 

TECHNICAL NOTES: Information less broad in scope but nevertheless of 
importance as a contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: Information receiving limited distribu- 
tion because of preliminary data, security classification, or other reasons. 

CONTRACTOR REPORTS: Scientific and technical information generated 

under a NASA contract or grant and considered an important contribution to 
existing knowledge. 

TECHNICAL TRANSLATIONS: Information published in a foreign 
language considered to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information derived from or of value to NASA 
activities. Publications include conference proceedings, monographs, data 
compilations, handbooks, sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION PUBLICATIONS: Information on tech- 
nology used by NASA that may be of particular interest in commercial and other 
non-aerospace applications. Publications include Tech Briefs, Technology 
Utilization Reports and Notes, and Technology Surveys. 


Details on the availability of these publications may fee obtained from: 


SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington D.C. 20546 



